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Abstract 

Ureolytically-driven calcium carbonate precipitation is the basis for a promising in-situ remediation method for 
sequestration of divalent radionuclide and trace metal ions. It has also been proposed for use in geotechnical 
engineering for soil strengthening applications. Monitoring the occurrence, spatial distribution, and temporal 
evolution of calcium carbonate precipitation in the subsurface is critical for evaluating the performance of this 
technology and for developing the predictive models needed for engineering application. In this study, we 
conducted laboratory column experiments using natural sediment and groundwater to evaluate the utility of 
geophysical (complex resistivity and seismic) sensing methods, dynamic synchrotron x-ray computed tomography 
(micro-CT), and reactive transport modeling for tracking ureolytically-driven calcium carbonate precipitation 
processes under site relevant conditions. Reactive transport modeling with TOUGHREACT successfully simulated the 
changes of the major chemical components during urea hydrolysis. Even at the relatively low level of urea 
hydrolysis observed in the experiments, the simulations predicted an enhanced calcium carbonate precipitation 
rate that was 3-4 times greater than the baseline level. Reactive transport modeling results, geophysical monitoring 
data and micro-CT imaging correlated well with reaction processes validated by geochemical data. In particular, 
increases in ionic strength of the pore fluid during urea hydrolysis predicted by geochemical modeling were 
successfully captured by electrical conductivity measurements and confirmed by geochemical data. The low level 
of urea hydrolysis and calcium carbonate precipitation suggested by the model and geochemical data was 
corroborated by minor changes in seismic P-wave velocity measurements and micro-CT imaging; the latter 
provided direct evidence of sparsely distributed calcium carbonate precipitation. Ion exchange processes promoted 
through NH 4 + production during urea hydrolysis were incorporated in the model and captured critical changes in 
the major metal species. The electrical phase increases were potentially due to ion exchange processes that 
modified charge structure at mineral/water interfaces. Our study revealed the potential of geophysical monitoring 
for geochemical changes during urea hydrolysis and the advantages of combining multiple approaches to 
understand complex biogeochemical processes in the subsurface. 



Background 

Remediation of subsurface radionuclides and trace metal 
contaminants is one of the U.S. Department of Energy's 
(DOE) greatest challenges for long-term stewardship [1]. 
Co-precipitation in calcium carbonate is one attractive in- 
situ remediation strategy for divalent radionuclide and 
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trace metal ions, such as 60 Co, 137 Cs, and 90 Sr [2-4]. The 
partitioning of these trace metals, e.g. 90 Sr, into the cal- 
cium carbonate lattice is partly a function of precipitation 
rate [5-7]. One mechanism to promote calcium carbonate 
precipitation is through urea hydrolysis which increases 
groundwater pH and alkalinity, resulting in greater 
mineral saturation with respect to calcium carbonate. 
Urea hydrolysis is catalyzed by the urease enzyme, 
expressed by many microorganisms in the subsurface in 
order to harvest nitrogen. Figure 1 illustrates the 
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Figure 1 Conceptual diagram of ureolytically-driven calcium carbonate precipitation approach for remediation of 90 Sr contaminated 
geologic media. A. Microbially catalyzed hydrolysis of urea. B. Cation exchange and calcium carbonate precipitation. C. Continued precipitation 
of calcium carbonate. Figure adapted from [8]. 



conceptual model of the proposed reaction and contami- 
nant sequestration processes [8]. 
The reactions involved in this process include: 

1. Enzymatically catalyzed urea hydrolysis produces 
NH 4 + and HC0 3 and raises pH: 

H 2 N(CO)NH 2 +H + + 2H 2 0 -> 2NH + 4 + HCO~ (1) 

2. NH 4 + promotes desorption of cations M 2+ (e.g. Ca 2+ 
and Sr 2 *) from grain surfaces: 

Solid - (A4 2+ ) + 2NH + 4 0 Solid - [NH + 4 ) 2 + M 2+ (2) 

3. HCO 3 promotes precipitation of calcium carbonate 
and co-precipitation of 90 Sr: 

x 90 Sr 2+ + (1 - x)Ca 2+ + 2HCO3 & Ca ( i_ x) 90 Sr x CO s + H 2 0 + C0 2 (3) 

The net reaction leads to the production of carbonate 
minerals containing Sr: 

Solid-{x c,0 Sr 2+ )+H 2 N{CO)NH 2 +{l-x)Ca 2+ +2H 2 0 — Ca (l - x) Sr x C0 3 +Solid-{NH + 4 ) 2x (4) 

Compared to direct injection of carbonate solution that 
could lead to clogging near the injection wells due to 
rapid precipitation of calcium carbonate, microbially 
facilitated calcium carbonate precipitation relies on the 
production of carbonate via urea hydrolysis, allowing 
urea to transport further into the subsurface before sig- 
nificant production of carbonate occurs. This could miti- 
gate the risk of wellbore clogging and allow treatment of 
a larger area while using a single injection point. In addi- 
tion, the NH 4 + produced during ureolysis can exchange 
with cations (including radionuclides such as 90 Sr 2+ ) pre- 
viously sorbed to the subsurface materials (Eq. 2), thereby 
promoting metal sequestration in a more stable form 
(precipitation instead of sorption) (Eq. 3). This remedia- 
tion approach is particularly attractive for environments 
where the groundwater is already saturated or oversatu- 
rated with respect to calcium carbonate, and thus the 



potential for releasing trapped metals through calcium 
carbonate dissolution over the long term is minimal. The 
arid western U.S. is an example of such an environment 
[9]. In addition to metal sequestration, enhanced calcium 
carbonate precipitation has many other potential applica- 
tions, such as providing cementation for soil strengthen- 
ing [10], and a long term sink for atmospheric carbon 
associated with C0 2 storage/sequestration within geolo- 
gical media [11]. 

Although studies have demonstrated the potential use- 
fulness of enhanced calcium carbonate precipitation for 
subsurface applications [10,12-15], one critical issue yet 
to be resolved is the monitoring of the calcium carbonate 
precipitation process in the subsurface to track its occur- 
rence, spatial distribution, and temporal evolution. Such 
monitoring is necessary in order to evaluate its perfor- 
mance. In recent years, near surface geophysical methods 
have been used to provide high resolution imaging of 
subsurface properties in a minimally invasive manner, 
and these methods hold great promise for our ability to 
monitor subsurface perturbations [16,17]. Due to its sen- 
sitivity to mineralogy and pore fluid chemistry, electrical 
methods, e.g. complex conductivity, are among the most 
explored. 

Complex conductivity (<7*) measures the frequency 
dependent electric conduction behavior of a porous 
medium. The measured o*(co) of a porous medium can 
be represented as, 

G*{cD) = ol{cD) + ioll{co). (5) 

where co is the angular frequency, & is the measured 
real part of O*(co), being the conduction (energy loss) 
component; a" is the measured imaginary part of 0*(a>), 
being the polarization (energy storage) component and 
i = V— 1 • At low frequencies (< 1000 Hz), electrical 
charge transport in saturated porous media is deter- 
mined by two processes: the electrolytic conductivity 
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(<7 e /), representing conduction via interconnected fluid- 
filled pore space (a purely real term), and a complex 
interfacial conductivity (<7* nt ), representing conduction 
and polarization processes occurring near the grain/elec- 
trolyte interfaces. The c e i is dependent on the conduc- 
tivity of the electrolyte (<J W ) saturating the porous media 
[18] whereas o* int represents an interfacial charge con- 
duction and polarization behavior within the electrical 
double layer (EDL) at the grain/electrolyte interface 
dominated by a diffusive mechanism associated with ion 
migration to and from the grain surfaces [19-21]. The 
EDL is composed of (1) an inner fixed (Stern) layer with 
more or less structured charges adsorbed on the mineral 
surface through complexation and (2) an outer diffuse 
layer with more randomly distributed charges held 
weakly through electrostatic forces. Various factors 
could affect the interfacial charge conduction and polar- 
ization behaviors including mineralogical composition, 
specific surface area, particle size distribution as well as 
EDL properties, such as charge density, mobility and 
thickness. The EDL properties are influenced strongly 
by mineralogy and pore fluid properties, such as ionic 
speciation, concentrations, pH, etc. Other less explored 
factors that could change interfacial EDL properties 
include ion exchange behavior as well as the presence of 
organics [22]. 

When considering the polarization behavior at the 
grain/fluid interfaces of natural sediments, the inner 
Stern layer of the EDL structure is normally viewed as a 
relatively stable plane with primarily tangential move- 
ment of charges along the interface [21]. Charge move- 
ment normal to the interface is more difficult due to the 
higher potential barrier, the so-called zeta potential, that 
charges have to overcome in order to move away from 
or closer to the mineral surface [21]. Compared to the 
charges in the fixed layer, charge movement in the dif- 
fuse layer can occur more easily in either tangential or 
normal directions under external potential field due to 
the much weaker electrostatic forces. The polarization 
magnitude at the interface is primarily determined by 
the EDL properties (e.g. charge density, mobility). 
Alterations in pore fluid chemistry modify surface com- 
plexation structure, resulting in the changes of the EDL 
properties, thus polarization behavior [23]. This model 
has been successfully used to explain the polarization 
behavior of natural sediments. For example, under the 
condition of fixed solid matrix properties (mineralogy, 
surface area/particle size), previous research has shown 
that imaginary impedance (equivalent to a") changes in 
response to variations in ion concentrations in the pore 
fluid in Berea sandstone samples [24]. Other researches 
also demonstrated the sensitivity of electrical signals to 
changes in pore fluid chemistry and/or precipitation of 
new mineral phases due to amendment injection and 



mineral precipitation at both laboratory and field scales 
[17,25-29]. For example, Williams et al. [25] performed 
a laboratory biostimulation experiment where time-lapse 
complex conductivity, seismic, and various geochemical 
measurements were collected. They showed that 
changes in complex conductivity and seismic amplitude 
correlated with the onset and spatial distribution of 
microbially-mediated iron and zinc sulfide precipitation. 
Slater et al, [30], Ntarlagiannis et al, [27] and Personna 
et al [31] also showed the sensitivity of complex conduc- 
tivity to metal sulfide precipitation/dissolution. These 
experiments demonstrated the potential of geophysical 
methods, particularly electrical measurement techniques, 
for the monitoring of mineral precipitation in porous 
media at the column scale. Recent studies have also illu- 
strated the potential of geophysical methods to track 
remediation processes at the field-scale [16,17,28]. 

Particularly relevant to this study, Wu et al. [29] 
showed that abiotically induced calcium carbonate preci- 
pitation produced detectable geophysical signals during a 
laboratory column experiment. In that study, extensive 
and rapid calcium carbonate precipitation was induced 
through mixing of calcium and carbonate solutions in a 
glass bead packed column. Complex conductivity mea- 
surements collected daily during the experiment revealed 
a significant increase in the polarization magnitude, 
which coincided with the accumulation of calcium carbo- 
nate precipitates within the column and demonstrated 
clearly the sensitivity of the electrical signal to calcium 
carbonate precipitation. The polarization response 
observed was interpreted to be due to a large enhance- 
ment of the surface charge of the solid phase due to sub- 
stantial calcium carbonate precipitation. However, 
because the experiment was conducted abiotically with a 
relatively high reaction rate and within an artificial med- 
ium, the sensitivity of geophysical methods to calcium 
carbonate precipitation at lower, field relevant rates and 
in the presence of field derived sediments warrants 
further evaluation. 

In this study, we conducted laboratory flow-through 
column experiments to explore sensitivity of geophysical 
methods to ureolytically-driven calcium carbonate preci- 
pitation under field relevant conditions with natural sedi- 
ments, groundwater and microbial consortium. Primary 
mechanisms by which microbially enhanced calcium car- 
bonate precipitation through urea hydrolysis could 
change both the conduction and polarization compo- 
nents of the complex conductivity through (1) changing 
the ionic strength/speciation and pH of the pore water 
(and thus the surface charge properties of the solid), (2) 
increasing the total surface area of the solid matrix 
through the precipitation of calcium carbonate, and (3) 
changing the hydraulic conductivity/permeability as a 
result of changes in the pore geometry (if significant 
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calcium carbonate precipitation occurs). In addition to 
contributions to the increases in ionic strength, NH 4 + 
produced during urea hydrolysis (Eq. 1) could engage in 
ion exchange with the sediments (Eq. 2), thereby provid- 
ing additional polarization effects through altering sur- 
face charge structure by (1) changing surface charge 
density through non-stoichiometric charge exchange and 
(2) altering complexation structure between surface sites 
and absorbed species. 

In addition to electrical signals, seismic wave propaga- 
tion velocities have proven to be moderately sensitive to 
bioprecipitate formation with previous studies document- 
ing changes in P-wave attenuation during FeS precipita- 
tion [25,32] and significant increases in S-wave velocity 
during CaC0 3 precipitation induced by microbial urea 
hydrolysis have also been observed [13]. Generally, 
changes in seismic velocities are due to stiffening of the 
granular framework, a process that should be captured 
by existing contact cement theories (CCT) [33]. Replace- 
ment of pore fluids by a solid mineral phase at weak 
grain contacts results in an increase in both bulk and 
shear modulus during the precipitation process. Second- 
ary changes in seismic properties related to poroelastic 
phenomenon, including so-called double porosity effects 
[34], may also be present in this class of systems due to 
heterogeneous precipitation or micro-porosity within the 
precipitate phase. 

In addition to geophysical monitoring, ID reactive 
transport modeling using TOUGHREACT ([35] and 
references therein) was carried out to simulate reactions 
occurring during the experiment. Geochemical processes 
considered in the simulations included ureolysis together 
with calcium carbonate precipitation, ion exchange, and 
nitrification. The effect of dissolution and precipitation of 
minerals other than calcium carbonate was assumed to 
be negligible during the relatively short time frame of the 
experiment. 

Dynamic synchrotron x-ray micro-CT measurements 
were conducted to characterize precipitate distribution to 
inform analysis of observed geophysical signatures. X-ray 
micro-CT is a non-destructive 3D imaging methodology 
with micron scale resolution which has been applied to 
the characterization of porous solids since the mid-1990s 
[36,37] . By dynamic, we refer to multiple CT images cap- 
tured over a period of time to allow process monitoring 
[38]; in this case, the process is the formation of CaC0 3 
precipitates. The micro-CT technique is directly analo- 
gous to medical CT scanning with the exception of the 
greatly enhanced resolution, smaller imaging target, and 
a different x-ray source. In the case of precipitate genera- 
tion within a porous medium, a solid mineral phase is 
replacing the pore fluid, a change easy to observe using 
micro-CT methods. In this particular study, we use 
dynamic micro-CT data qualitatively to characterize the 



spatial pattern of precipitation on the pore scale. The key 
questions that we address include the spatial distribution, 
morphology, and amount of precipitates formed. 

We anticipated coupling between and uncertainty asso- 
ciated with several processes, including: the magnitude 
and spatial distribution of calcium carbonate precipitates; 
the clogging of pore throats and concomitant change in 
flow properties associated with precipitation [39]; and 
temporal variation of geophysical signatures associated 
with different phases of the reaction processes. The com- 
bination of geophysical datasets, reactive transport model- 
ing, micro-CT imaging and jointly interpreted monitoring 
and modeling datasets helped reducing ambiguities asso- 
ciated with these processes. Specifically, dynamic synchro- 
tron x-ray micro-CT measurements provided direct 
information about the magnitude and nature of calcium 
carbonate precipitation at the pore scale including its geo- 
metric arrangements within pores and spatial distribution 
throughout the column. Complex conductivity and seismic 
P-wave velocity measurements were used to gain under- 
standing about effective changes in geochemical and geo- 
mechanical conditions, respectively, averaged over the 
length of the column. Reactive transport modeling helped 
to understand the rates associated with the reaction 
network described in Equations 1-4 and to predict the 
volume of produced calcium carbonate. Together, the 
methods were used to gain an improved predictive under- 
standing of ureolytically-driven calcium carbonate precipi- 
tation in natural media and particularly, of the potential of 
geophysical methods for monitoring the relevant biogeo- 
chemical changes over space and time. 

Methods 

Column setup 

Two columns were fabricated for the experiment from 
polycarbonate plastic: a large multi-port column for elec- 
trical and seismic measurements and a second smaller col- 
umn for synchrotron micro-CT imaging. The larger 
column has a length of 20 cm and 7.6 cm inner diameter 
(I.D.) while the micro-column has a length of 12.7 mm 
and an I.D. of 6.35 mm. 

Alluvial sediments used in the experiment were col- 
lected from a background site north of Idaho Nuclear 
Technology and Engineering Center (INTEC) of the 
Idaho National Laboratory (INL). 90 Sr contamination was 
found in the subsurface below the INTEC facility. 
Groundwater used in the experiment was from the 
Vadose Zone Research Park (VZRP) at the INL. The 
VZRP subsurface has geological and hydrological proper- 
ties similar to those at the nearby INTEC. The VZRP 
groundwater is not radiologically contaminated but has 
naturally occurring nonradioactive Strontium that can be 
used to study the co-precipitation behavior during ureo- 
lytically-driven calcium carbonate precipitation. 
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For the large electrical/seismic column, the sediments 
were used after large pebbles (> 1 cm) had been removed. 
This sediment had been archived and dried for over ten 
years. Significant clay content ranging from 8% to 28% 
has been measured for core retrieved from the site, 
depending on locations and depth; illite and smectite are 
among the major clay minerals identified [40,41]. The 
dry material was saturated in VZRP groundwater (fluid 
conductivity, o w ~ 0.062 S/m) overnight and then wet 
packed into the column. A 6 cm section of coarse silica 
sand (particle size ~1 mm) was packed at the bottom of 
the column to facilitate homogeneous fluid flow into the 
sediment. In an effort to maximize homogeneity, the 
sediments were packed ~ 1 cm at a time and tapped 
from outside to allow settlement and compaction [42]. 
However, seismic data revealed a certain level of hetero- 
geneity along the column length (discussed below). The 
smaller micro-CT column was packed using the same 
protocol except that the sediment was sieved to include 
only particles less than 1 mm in size; this modification 
was chosen due to the small I.D. of the column (6.35 
mm). An initial hydraulic conductivity (K) measurement 
was carried out on the electrical/seismic column using 
the constant head method and a K value at 1.2E-5 m/s 
was measured. This value is within the range of values 
previously measured on core samples retrieved from 
boreholes drilled at the INL (varying from 1E-5 to 1E-9 
m/s depending on depth and location) [40]. 

The experiment was carried out in three phases: Phase 
I: groundwater was pumped through the columns to 
achieve an initial steady state as defined by stability of 
effluent conductivity; Phase II: groundwater amended 
with molasses (5 mg/L) was injected into the columns for 
~3 days to stimulate microbial growth and activity [15]; 
and Phase III: groundwater amended with urea (10 mM) 
was delivered into the columns to induce urea hydrolysis 
and calcium carbonate precipitation. The molasses con- 
centration was purposefully kept low to stimulate micro- 
bial activities at an appropriate level, (i.e. to prevent pore 
clogging by excessive precipitation or the transition into 
a sulfate reduction phase). 

For the large column experiments, complex conductiv- 
ity measurements were acquired using a four electrode 
configuration, with one electrode installed inside each 
end cap for current injection and two installed in the 
middle of the columns for measurements (Figure 2). Ag/ 
AgCl electrodes were used for both current injection and 
potential measurements. Spiral current electrodes inside 
the end caps were in direct contact with column material 
while potential electrodes were electrolytically coupled 
with the column through small plastic tubes installed on 
the column wall and filled with the groundwater saturat- 
ing the column. The electrodes recovered after the 
experiment had no visible alterations (e.g., mineral 



precipitation or dissolution of AgCl coating). For the 
measurements of the seismic signals, four pairs of P-wave 
seismic transducers (1 MHz) were installed along the col- 
umn length in an orientation 90 degrees from the Ag/ 
AgCl potential electrodes. Each transducer was coupled 
to the sediment through a 0.5 cm aluminum buffer rod 
which directly contacted the material within the column. 

A constant flow rate of 0.21 ml/min (-0.7 pore 
volumes per day with a measured porosity of 45%) was 
used throughout the experiment for the large column, 
which resulted in a residence time of -32 hours. A frac- 
tion collector (Foxy 200, ISCO) was used to collect efflu- 
ent samples on regular basis from the large column. 
Hydraulic conductivity measurements were conducted at 
the end of the experiment to determine the change in 
hydraulic conductivity from calcium carbonate precipita- 
tion. The injection rate for the micro-column was 
selected to match the fluid residence time (0.7 pore 
volumes/day) in the larger column and the same influent 
supply was used for both columns. 

Geophysical methods 

Electrical measurements were collected with a National 
Instruments (NI) dynamic signal analyzer (DSA, NI 
4461) using electrodes placed along the length of the 
large column (described above). A preamplifier was used 
to boost the input impedance to 10 9 Ohm to avoid cur- 
rent leakage into the measurement circuitry. Water col- 
umn calibration and repeatability tests indicate that 
errors were less than 0.3 mrads for the phase and 0.5% 
for conductivity when low frequency was used (< 500 
Hz). Each measurement is composed of a phase shift (0) 
and a magnitude (|o"|) component recorded relative to a 
precision reference resistor for forty frequencies spaced 
at equal logarithmic intervals from 0.1 to 1000 Hz. The 
real and imaginary parts of the complex conductivity 
represent the magnitude of the conduction and polariza- 
tion of the sample respectively, and were determined 
from the following equations: 

o' = \cr\ X COS0 (6) 

o" = \g\ x sin0 (7) 

In this study, we also measured P-wave velocity and 
amplitude using the ultrasonic transmission method at 
~1 MHz. At several points along the length of the col- 
umn, piezoelectric compression transducers mounted 
opposite each other allowed a through-column measure- 
ment of P-wave traveltime at regular time intervals dur- 
ing the experiment. The transducers were excited using 
a broadband ultrasonic pulser (Panametrics Model 
5075) and the transmitted signal was recorded via a 
digital oscilloscope (Tektronix Model TDS210). The 
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Figure 2 Schematics of the flow through column design. Black dots represent locations of the seismic transducers and electrodes are 
indicated in red. 



resulting transmitted waveforms were picked using a 
customized Graphic User Interface (GUI) and then cor- 
rected for zero-time and column shell delays. The times 
were then converted to P-wave velocities using the 
known transmission distance. 

Geochemical monitoring 

Geochemical analysis of effluent samples from the large 
column was conducted on regular basis. Electrolyte con- 
ductivity (a w ) and pH were measured near real time with 
a pH/ conductivity meter (Orion 5 star, Thermo Scientific). 
Effluent dissolved oxygen (DO) was measured with a YSI 
Model 5000 DO meter and was used as an proxy indicator 
of microbial activity during molasses injection and was 
measured periodically to ensure sustained microbial activ- 
ity afterwards during urea injection. Major anions (CL, 
S0 4 2 " and N0 3 ") were measured with ion chromatography 
(IC, Dionex) and major cations were measured with 
Inductively Coupled Plasma Mass Spectrometry (ICP-MS, 
Perkin Elmer). Alkalinity was measured using an acid 



titrimetric method (Mettler Toledo DL 50 Graphix titra- 
tor). Effluent ammonia and urea concentration were also 
quantified using IC analysis. Due to the automated sam- 
pling procedure with the fraction collector and the equili- 
brium of the sample fluid with atmospheric pC0 2 , partial 
NH 4 + oxidation (production of N0 3 ") in the effluent and 
premature calcium carbonate precipitation in the influent 
container were observed and both issues were accounted 
for in the modeling studies. 

Micro-CT imaging 

A 3 rd generation synchrotron, the Advanced Light Source 
(ALS) at Lawrence Berkeley National Laboratory was 
used as the x-ray source for micro-CT imaging. Synchro- 
tron x-ray sources are attractive for micro-CT studies 
because of their extremely high flux and parallel beam; 
the former feature allows for use of a monochromator 
avoiding beam hardening artifacts. All imaging for our 
experimental effort was conducted at ALS beamline 
8.3.2, which is positioned on a super-bend magnet and 
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provides access to x-ray energies between 8 and 41 keV. 
The resulting x-ray attenuation image stack is sensitive 
to changes in localized density or atomic number, thus 
allowing discrimination between solid, fluid, and gas 
phases as well as some differences in solid mineral 
phases. 

The micro-CT acquisition protocol includes careful 
mounting of the experimental subject on a rotating stage 
and acquisition of a sequence of x-ray projections over a 
range of 180 degrees. X-rays passing through the sample 
are converted to visible light at a scintillator, and the 
resulting light passes through a traditional microscope 
optics chain and is recorded using a high-resolution CCD 
(Cooke PCO4000). This sequence of radiographs is then 
inverted using a filtered back-projection algorithm (Octo- 
pus, InCT Inc.) to generate a 3D volume of x-ray absorp- 
tivity. A 2x optical objective (Mitutoyo) was used during 
scan yielding final 3D reconstructed volumes with cubic 
voxels of -4.47 microns. The entire length of the column 
was scanned for both the baseline and repeat imaging 
runs. 

Destructive sampling 

The columns were destructively sampled after the experi- 
ments for multiple characterization purposes. The large 
geophysical column was sub-cored with polycarbonate 
tubing for micro-CT imaging and comparison with the 
miniature column to characterize geometry and spatial 
distribution of the precipitates. Both pore fluids and fine 
sediment samples from the column were characterized for 
cell density after staining with acridine orange and direct 
counting using a Zeiss Axioskop fluorescence microscope. 
In addition, cell adherence to sediment particles was deter- 
mined through washing with CaCl 2 and Tween 80/pyro- 
phosphate solutions to neutralize cell surface charge and 
promote cell detachment [43] and cells in supernatant 
wash fluids were counted after staining. Recent research 
found that ureolytic activity on a per gram or per ml basis 
was much higher on the sediments than in the fluid [8]. 
Cell detachment analysis provides information of relative 
abundance of cells attached to sediments comparing to 
those in the fluid, and can be used to evaluate relative 
ureolytic activity of these cells. 

Reactive transport modeling 
Model Setup and Approach 

The numerical code TOUGHREACT ([35] and references 
therein) was used to simulate reactions occurring during 
the flow through column experiment. A 15-day period 
was simulated, starting at the time of urea-amended 
groundwater injection (Phase III). The initial sediment 
pore water was taken as the measured water composition 
at the column outlet after the initial equilibration period 
(Phase I), but prior to urea injection (Table 1). However, 



the elevated measured pH (-8.6) at the column outlet at 
the end of Phase I indicated partial equilibration with 
atmospheric C0 2 upon sampling and supersaturation of 
the solution with respect to calcium carbonate. To better 
represent conditions within the column, the initial total 
carbonate concentration and pH within the column were 
set to reflect equilibrium with calcium carbonate, which 
yielded pH and total dissolved carbonate values (Table 1) 
more in line with pH values determined for VZRP 
groundwater. 

The pH measured in the influent reservoir was higher 
(~ 8.2) than that observed in the VZRP (-7.6) where the 
groundwater used in the experiment was acquired, also 
indicating some re-equilibration with atmospheric C0 2 
in the reservoir. In this case, the composition of the 
injected water (column inlet) was specified using the 
actual measured pH and species concentrations in 
the influent reservoir, but correcting the total dissolved 
carbonate concentration to yield ion charge balance 
(Table 1). 

Geochemical processes considered in the simulations 
included ureolysis together with calcium carbonate preci- 
pitation, ion exchange, and nitrification, as further dis- 
cussed below. The effect of dissolution and precipitation 
of minerals other than calcium carbonate (such as silicates 
and aluminosilicates) was assumed to be negligible during 
the relatively short time frame of the experiment on the 
basis that the reaction rates of these minerals are much 
slower than calcite reaction rates [44]. Calcium carbonate 
was allowed to precipitate in the sand pack, but this part 
of the column was otherwise assumed unreactive, with 
essentially no ion exchange capacity. 

The model domain was set as a ID cylindrical column 
with dimensions and hydraulic properties described 
above. A fixed influx rate (0.21 mL/min) was applied at 
the modeled column inlet. The column was discretized 
into 205 grid blocks at regularly spaced intervals of 1 
mm. A sequential-iterative (transport/reaction) method 
was implemented to reproduce accurate concentration 
time profiles, using a maximum time step of 416 s (-0.8 
x Courant). Tests using finer time and space intervals 
indicated that this set of parameters provided a reason- 
able compromise between simulation speed and accu- 
racy. To allow for partial re-equilibration of effluent 
water with atmospheric C0 2 , the modeled column outlet 
was set at a constant pC0 2 of 10~ 3 ' 3 bar, a value estimated 
(by trial and error) to yield best results; C0 2 was not 
allowed to diffuse back into the column. Automatic cali- 
bration of other model parameters was performed using 
PEST [45] on a Unix cluster, including the total enzyme 
concentration driving the hydrolysis of urea (including 
EH species in Eq. 8, 9 and 10 shown below), a similar 
concentration term driving the rate of nitrification (Q, i0 / 
Y bio in Eq. 13 shown below), the calcium carbonate 
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Table 1 Water compositions used in simulations 
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IllllUcill LLM 1 IJJUbl ULM 1 UUIIIIL) [Jlldbcr III 
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mM 
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3.03 


ivig 


mM 
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0.487 


A 1 

Al 


mM 


O.zo X 1 0 


0. 1 5 X 1 0 


Si 


mM 
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k" 
t\ 


mM 


n 09.1 

U.Uo/ 


n 

U.UDD 


Ca 


mM 


1.05 


1.00 


Sr 


mM 


2.82 x 10~ 3 


2.83 x 10~ 3 


CI 


mM 


2.63 


2.63 


N0 3 


mM 


44.5 x 10~ 3 


18.1 x 10~ 3 


S0 4 


mM 


0.213 


0.206 


PH 




7.62 (a) 


8.23 


EC0 3 (a) 


mM 


3.23 (a) 


4.21 (a) 


PC02 


bar 


10 -2.37 (a) 


1Q-3.01 (a) 


Urea 


mM 


nil (10" 1 °) 


10 


a Calculated (see text) 



precipitation rate constant, and the ion exchange capacity 
of the sediments. 
Ureolysis kinetic rate law 

Ureolysis was modeled as an enzymatic reaction, with a 
rate law taking into account the effect of pH on enzyme 
(urease) activity as well as product inhibition by NH 4 + [46]: 



R = k- 



[EH][Urea] 



(JC M + [Urea 



[NHJ] 



10-P H K 2 



(8) 



10-P H / 



In this equation, brackets indicate concentrations, EH 
represents the enzyme driving the hydrolysis of urea, K M 
and K P are the half-saturation and inhibition constants, 
respectively, and k is a rate constant. Kj and K 2 are disso- 
ciation constants representing the enzyme protonation 
and de-protonation reactions, respectively. This rate law 
was implemented into TOUGHREACT, as described pre- 
viously [47], using a standard Michelis-Menten rate law 
with an inhibition term: 



R = fe[EH] 



[Urea] 



K M + [Urea] K P + [NH+] 



(9) 



A primary species, EH, was added to the input ther- 
modynamic database to represent the enzyme, together 
with two additional derived species, EH 2 + and E", 
expressing the protonated and deprotonated enzyme 
species: 



EH+ EH + H + log(Ki) 25 o C = -6.121 



E" + H + ±* EH log(K 2 ) 2 5°c = 7-896 



(10) 



(11) 



It can be shown that this approach is a mathematical 
equivalent to using Equation 8. Values of K 2 and K 2 
were taken from Fidaleo and Lavecchia [46]. The total 



concentration [EH] + [EH 2 + ] + [E~] was assumed con- 
stant and calibrated (10~ 10 mol urea se/L)- Other rate law 
parameters were taken from Fidaleo and Lavecchia [46] 
as follows (at 25°C, with k re-expressed using molar 
concentrations): K P = 1.22 x 10 2 mol/L, K M = 3.21 x 
10~ 3 mol/L, and k = 146.4 mol urea mol urease _1 s" 1 . 
Calcium carbonate precipitation 

In the model, calcium carbonate (CaC0 3 ) was allowed to 
form as an ideal solid solution with strontianite (SrC0 3 ) 
(i.e., assuming equality of activity and mole fraction). The 
solid solution was set to precipitate under kinetic con- 
straint using a rate constant calibrated from the experi- 
mental data. A standard rate law derived from Transition 
State Theory (TST) was implemented [48], in which the 
rate constant is multiplied by (1-Q/K), with Q and K 
representing the ion activity product and equilibrium 
constant, respectively, of the calcium carbonate/strontia- 
nite solid-solution precipitation reaction. Nucleation was 
not specifically modeled, and the rate constant was cali- 
brated as a bulk average rate (4.2 x 10" 10 mol/s) incorpor- 
ating surface area effects. Although this approach is quite 
approximate, it yielded satisfactory results. More elabo- 
rate, mechanistic precipitation models were deemed 
unwarranted given that the uncertainty of the calcium 
carbonate precipitation rate is large and partly carried 
through the calibration of other parameters. 
Ammonium oxidation 

The oxidation of NH 4 + was considered in the simula- 
tions to account for an increase (-0.02 mM) in nitrate 
concentrations observed at the onset of ureolysis. The 
oxidation was attributed to nitrification, which was 
approximated by the following overall reaction (assum- 
ing no NO 2 buildup): 



NH; + 20 



2(aq) 



NO3 +H 2 0 + 2H + 



(12) 
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The overall rate was described using Monod kinetics: 

Mmax r C NH + 4 C Q 2 H) , . 

KnHX = ^7— C bio- — " — (13) 

mo &NH + 4 + <-nh+ K o 2[ac1) + <-o 2M 

Half-saturation constants K NH4 + (1.48 x 10" 5 mol/L) 
and Ko2(aq) (2.41 x 10~ 5 mol/L), and the rate (i max (9.53 x 
10" 6 s" 1 ) were taken from the literature [49]. The Cb io IY bio 
term was assumed constant and was calibrated (10~ 15 
mol/L). The dissolved oxygen concentration at the col- 
umn inlet (C0 2 (aq) in Eq. 13) was set to 2 x 10-3 mol/L 
(~ 6 mg/L, close to atmospheric), in the range of mea- 
sured concentrations. These concentrations dropped by 
about half in the column effluent. Note that the experi- 
ment was carried out specifically under oxic conditions 
to avoid potential complications arising from oxygen 
depletion by microbial reactions. Because oxygen-limited 
conditions were not observed, and amounts of produced 
nitrate were consistent with the amount of oxygen 
decrease, other microbially-driven oxygen-consuming 
reactions (e.g. organic matter oxidation) did not appear 
to be limiting and thus were not included into the model 
for simplicity. 
Ion exchange 

Ion exchange was implemented for Na + , K + , NH 4 + , Ca 2+ , 
Mg 2+ , and Sr 2+ , using a generic exchanger (i.e., no speci- 
fic minerals are associated with exchange). Ion selectivity 
coefficients for INL sediments from Baker et al. [50] were 
used, except for K + , for which no data were reported by 
these authors, and for which the same selectivity as for 
Na + was assumed. The implemented exchange model fol- 
lowed the Gapon convention [51] with Ca 2+ as the refer- 
ence cation. The sediment cation exchange capacity was 
calibrated, yielding a value (11.9 cmol/kg) in line with 
that measured by Baker et al. [50] (21 cmol/kg) and other 
data cited by these authors (0.9-45 cmol/kg) for INL site 
materials. 

Results 

Aqueous geochemistry and microbial analysis 

Effluent geochemical data over the three phases of the 
experiment are shown in Figure 3. 

Major ions monitored during phase I include Ca 2+ 
(Figure 3A); Sr 2+ (Figure 3B); and S0 4 2 ~ (Figure 3C). S0 4 2 ~ 
was monitored as an indicator of whether the system pro- 
ceeded into sulfate reduction. During phase I, when site 
groundwater was pumped into the column, the major ions 
of concern (Ca 2+ , Sr 2+ and S0 4 2 ~) started at higher values 
compared to those measured later in this phase, when 
concentrations essentially matched the column inlet water. 
This was attributed to the effect of salt dissolution and ion 
desorption from the column matrix material. Alkalinity 
was relatively stable throughout this phase (Figure 3E). 
After ~4 days of equilibration with site groundwater, the 



column was stabilized with respect to most of the ions 
of interest including Ca 2+ , Sr 2+ , alkalinity and S0 4 2 ~ 
(Figure 3). 

During phase II, groundwater containing 5 mg/L 
molasses was injected into the column to stimulate 
microbial activity. DO was monitored on a daily basis 
during molasses injection and was used as an indicator 
of microbial activity. The effluent DO was measured at 
3.5, 2.2 and 1.7 mg/L during the first three days follow- 
ing molasses injection. Comparing to the constant influ- 
ent value at ~6 mg/L, the decrease of DO indicated the 
increased microbial activity. Phase III commenced when 
the effluent DO level dropped to ~ 1 mg/L to prevent 
the system from entering denitrification and sulfate 
reduction conditions. No other major geochemical 
changes were observed during this phase. 

For the final experimental Phase III, 10 mM urea 
amended groundwater (no molasses) was introduced into 
the column to initiate urea hydrolysis and calcium carbo- 
nate precipitation. Effluent geochemical monitoring 
revealed a decrease in urea concentration (compared to 
the influent level, Figure 3D), an increase of alkalinity 
(Figure 3E) and the production of NH 4 + (Figure 3F), all 
indicative of the occurrence of urea hydrolysis. However, 
the relatively high effluent urea concentration (average 
-8.7 mM compared to influent at 10 mM, Figure 3D) 
indicated that only -13% of the urea added was hydro- 
lyzed. This hydrolysis rate was stable during phase III. A 
peak concentration was observed at ~ day 16 for Ca 2+ 
and Sr 2+ at the beginning of phase III, coinciding with 
the production of NH 4 + . This is likely due to ion 
exchange of the NH 4 + for sorbed Ca 2+ and Sr 3+ at 
mineral surfaces (Eq. 2). Following the short duration 
peaks, the concentrations of Ca 2+ and Sr 2+ decreased gra- 
dually to steady levels, with a 28% decrease for Ca 2+ 
(from 1.08 mM to 0.78 mM) and 20% decrease for Sr 2+ 
(from 3 \iM to 2.4 (iM) compared to influent concentra- 
tions during Phase III. The trend for alkalinity was simi- 
lar to Ca 2+ and Sr 2+ with a short period of peak value 
upon the start of urea hydrolysis followed by a gradual 
decrease from 4.4 meq/L to 3.7 meq/L (a 16% decrease). 
The stoichiometric ratio between decreased alkalinity and 
Ca 2+ is ~ 2:1, consistent with precipitation of calcium 
carbonate according to Equation 3. Although an increase 
of pH is expected from urea hydrolysis according to 
equation 1, the pH of the effluent was stable at ~ 8.5 
throughout the experiment (data not shown), partially 
due to pH buffering by the sediments. As mentioned pre- 
viously, this pH value at 8.5 was higher than those 
observed in the VZRP where the groundwater used in 
this experiment was acquired, probably due to C0 2 out 
gassing to the atmosphere. S0 4 2 " concentrations were 
stable throughout Phase III, indicating no occurrence of 
sulfate reduction (Figure 3C). 
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Fluid conductivity (a w ) is sensitive to changes in total 
dissolved solids (TDS) and a distinctive pattern of o w 
change was observed during the different phases of the 
experiments that correlated well with geochemical data 
(Figure 3). The initial effluent o w (~ 0.07 S/m) from the 
column was higher than that of the groundwater (0.062 SI 
m) due to the accumulation of dissolved and desorbed 
ions in the pore water before reaching steady state. o w 
reached steady state in ~ 4 days and stayed relatively stable 
throughout the remainder of Phase I (Figure 3). During 
Phase II, the small concentration of molasses (5 mg/1) 
added had minimal impact on <J W . During phase III (urea 
injection) the fluid conductivity initially increased rela- 
tively rapidly from 0.064 S/cm to 0.0714 S/m (-12% 
increase), and then slowly decreased although it still 
remained elevated at the end of this phase. 

The total volume and surface area of calcium carbonate 
precipitated was estimated. A total of 4.8 liter of urea 
amended solution was pumped through the large column 
during Phase III. Based on the Ca 2+ concentration 



decrease in the effluent, the amount of calcium carbonate 
precipitated was estimated at an average of 0.12 milli- 
mole per liter of solution, resulting in a total of 0.058 g 
CaC0 3 precipitation during the experiment. This 
resulted in a total volume of 0.021 cm 3 based on a den- 
sity of 2.7 g/cm 3 for calcium carbonate. Assuming a uni- 
form distribution of the precipitates within the site 
sediments, this is equivalent to a negligible volume (~ 
0.007%) of the total pore space and one would also infer 
a negligible contribution to the change of total surface 
area. Hydraulic conductivity measurements made at the 
end of the experiment showed no changes relative to 
initial conditions, supporting the interpretation of limited 
calcium carbonate precipitation. 

Microbial cell densities in sediments (fine portion) and 
pore fluid samples were counted under a fluorescent 
microscope after staining with acridine orange and fixa- 
tion. In a dilute clay suspension (~ 0.1 g/1), the cell den- 
sity in the sediment is much lower (1.3E5) than that of 
the pore water (1.9E6). Cell adherence analysis showed 
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that most of the cells that appeared attached to the sedi- 
ments detached easily. This indicates either poor affinity 
of the microbial cells to sediment particles or the possibi- 
lity that these cells were from the fluid but artificially 
attached/appeared attached to sediments during sample 
preparation/fixation. 

Electrical Geophysical response 

Complex conductivity measurements (both & and a") 
acquired during the experiment revealed electrical signal 
changes associated with the physicochemical transfor- 
mations during enzymatically catalyzed urea hydrolysis 
(Figure 4). 

Figure 4 shows the changes of & and a" during the 
experiment calculated based on the phase shift (&) and 
magnitude (|cr|) measurements, as shown in equations 6 & 
7. As previously discussed, & is controlled primarily by 
charge transfer through interconnected fluid filled pore 
spaces in natural sediments. Our data shows that the 
changes of & (at 1 Hz shown in Figure 4A) followed the 
changes of o w throughout the experiment, indicating the 
dominant control of o w on bulk conductivity. During 
Phase I, when groundwater was pumped into the column 
for initial equilibration, both & and a w experienced sharp 
decreases before reaching a steady state. No significant 
changes were observed during molasses injection (Phase 
II). Both & and o w increased concurrently following urea 
injection in Phase III. After the first few days in Phase III, 
however, both & and o w began gradually decreasing dur- 
ing the course of urea injection, although the values 
remained elevated relative to the end of Phase I and to 
Phase II. 

As discussed previously, o" represents charge polariza- 
tion associated with the EDL and is affected by the 
charge structure of the specific surfaces of the solid 
phases determined by mineralogy as well as pore fluid 
properties. In our experiment a" increased step-wise 
throughout the experiment at all frequencies (Figure 4B). 
Percent changes of s" are plotted in Figure 4B for easier 
comparison between the three frequencies (1, 10 and 
100 Hz). During Phase I, a" increased steadily for ~ 7 
days before it leveled out for the reminder of Phase I and 
II. Although changes of electrical signatures associated 
with the introduction of microbial cells have been 
observed [52], no clear evidence of microbial effect was 
observed during phase II when the microbial activity was 
stimulated with molasses. Additional increase of a" 
occurred at the beginning of Phase III in parallel with the 
increase of o w due to urea hydrolysis. Once urea hydroly- 
sis and a w reached steady states, a" was stabilized and 
was maintained at a relatively constant level throughout 
the remainder of Phase III (Figure 4B). The actual values 
of the phase and imaginary conductivities before and 
after the experiments are listed in table 2. A comparison 



of o" trend at three different frequencies also shows a 
more significant increase at low frequencies (30 ± 7%, 
18 ± 4% and 7 ± 3% at 1, 10 and 100 Hz, respectively). 

Seismic results 

P-wave seismic transmission measurements acquired 
during the column experiment indicated only small 
changes in seismic properties and no direct signature of 
the stimulation of urea hydrolysis and calcium carbonate 
precipitation. Figure 5 shows the results of seismic mea- 
surements at the 4 discrete depth levels. Levels 1, 2, and 
3 are measuring property changes in the native sediment 
while level 4 sampled the sand at the bottom of the col- 
umn. Baseline velocities varied between 1600 and 1450 
m/s with variations attributed to packing heterogeneity 
and gravel distribution. As can be seen in panel A of Fig- 
ure 5, P-wave velocities displayed an initial gradual 
reduction for the transducer levels sampling the INL 
sediment during Phase I, while the sand showed an 
almost constant velocity of 1590 m/s. The maximum 
velocity reduction is ~ 60 m/s (~ 4%) at level 2 (Panel B). 
No signals correlated with either the molasses injection 
(Phase II) or the urea injection (Phase III). 

Synchrotron micro-CT Imaging Results 

Synchrotron micro-CT imaging was conducted on the 
small column after the initial packing and then again after 
2 weeks of urea injection. While a small amount of soil 
disturbance was noted at several points in the column, the 
granular pack was largely undisturbed. Positioning repeat- 
ability on the order of 20 microns was observed and man- 
ual registration was utilized to translate the volume for 
comparison. 

Figure 6 shows micro-CT imagery from 2 regions within 
the column before and after urea treatment. Several loca- 
lized regions with high density precipitates are visible in 
the repeat scan; these patches are typically 200-500 
microns in width, encompassing several grains. Within 
these regions, the precipitates manifested as grain coatings 
with no preference shown for grain-to-grain contacts. 
However, the number of precipitate patches is relatively 
small; the entire column contained only 15 zones with 
mineralization resolvable using micro-CT. No patches 
with similar characteristics were visible in the baseline 
images. There was no apparent preferential distribution of 
precipitation patches within the column overall, i.e. they 
occurred throughout the column rather than clustering 
near the inlet or column boundaries. 

Due to uncertainty as to whether the micro-column 
experiment was capturing the same precipitation pro- 
cesses as the larger column, a sub-core was extracted 
from the large column after completion of the urea treat- 
ment by insertion of rigid polycarbonate tubing. This 
sub-core was scanned at the same resolution as the 
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Figure 4 Complex conductivity measurements (A) & at 1 Hz and (B) c" at 1, 10 and 100 Hz during the experiment, both plotted 
together with the fluid conductivity data (black dots). 



micro-column (4.47 micron voxels) for direct comparison 
of precipitate architecture. Figure 7, panel A shows a CT 
cross-section of the entire sub-core; the large gaps in the 
sample are caused by the presence of small cobbles in the 
sediment which disrupted the fabric during coring, likely 
destroying some evidence of mineralization. However, 
careful examination of less disrupted zones within the 



core revealed high density features (Figure 7 B) similar to 
the observation of the small column (Figure 6 B, D). 

An attempt was made to directly examine a precipitate 
patch by embedding the core in epoxy and thin-section- 
ing for micro-XRD analysis; unfortunately the patches 
identified via CT could not be recovered due to sedi- 
ment fabric disruption during the embedding process. 



Table 2 Changes of the phase and imaginary conductivity during the experiment 



Frequency (Hz) 


&' (S/m) 




Percent change 


0 


(mRads) 


Percent change 


Before 


After 




Before 


After 




1 


5.1E-05 


6.7E-05 


30.6 


4.1 


5.4 


30.6 


10 


8.5E-05 


1 .0E-04 


18.2 


6.9 


8.1 


18.0 


100 


1 .2E-04 


1.3E-04 


7.7 


9.5 


10.3 


7.4 
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Reactive transport modeling results 

Several reactive transport simulations were run to isolate 
effects of ion exchange and ureolysis. The model was run 
suppressing ion exchange in one case, and suppressing 
ureolysis in another case. The measured trends of key 
constituents were well captured by the full model (Figure 
8), including sharp increases of NH 4 + from ureolysis (Eq. 
1) and of N0 3 " from NH 4 + oxidation (Eq. 12). Note that 
the sharp NH 4 + and N0 3 " concentration fronts could 
only be modeled by iterating between transport and reac- 
tion computations. Early effects of ion exchange were 
reproduced reasonably well for Ca 2+ and Sr 2+ , showing 
small initial peaks in concentrations followed by a steady 



decline probably due to calcium carbonate precipitation 
reducing the TDS level in pore fluid (Figure 8). As would 
be expected, the case without ion exchange shows a shar- 
per front for NH 4 + , no peaks in Ca 2+ and Sr 2+ concentra- 
tions, and much lower long term Sr 2+ concentrations in 
the effluent as less Sr 2+ is displaced by NH 4 + . The pH, 
which might have been expected to increase due to ureo- 
lysis (Eq. 1), was buffered by calcium carbonate precipita- 
tion (Eq. 3), as well as the imposed equilibration of the 
effluent with C0 2 at near atmospheric conditions. The 
amount of calcium carbonate forming in the column was 
predicted to be quite small. Note that the original inlet 
solution (groundwater) was predicted to be 
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Figure 6 Time-lapse synchrotron micro-CT images of precipitate formation. Panels A and C show two sections of the micro column before 
urea treatment. Panels B and D show the same two sections after 2 weeks of urea injection. Localized regions of grain-coating precipitation are 
visible. 



supersaturated with respect to calcium carbonate by 
about half a saturation index unit. As a result, precipita- 
tion is predicted to occur throughout the entire column, 
in part from ureolysis within the column and also 
because of supersaturation prior to contact with the col- 
umn sediments. The effect of ureolysis alone is apparent 
by comparing the results of simulations with and without 
ureolysis (Figure 9A and 9B), which show that 3 to 4 
times more calcium carbonate is predicted to precipitate 



when ureolysis is active, compared to the case when it is 
inactive. Note that the computed calcium carbonate 
amounts in the effluent (Figure 8) after 15 days are nearly 
double those predicted within the column (Figure 9) 
because the imposed pC0 2 at the column outlet (~10~ 3 ' 3 
bar, to reflect near-atmospheric conditions) is lower than 
within the column (~10~ 3 bar), resulting in additional cal- 
cium carbonate precipitation (e.g. 2HC0 3 " + Ca 2+ — > 
CaC0 3( s) + C0 2(g) t + H 2 0). 
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Figure 7 Selected micro-CT slice of sub-core extracted from large column. Panel A shows the cross-section of the entire sub-core while 
panel B shows a zoom of the small inset with a visible zone of localized mineralization. 



Discussion 

The production of NH 4 + and loss of urea showed that we 
successfully stimulated urease activity in our experiment 
following the provision of molasses for a short period of 
time. Although no other nutrients were supplied during 
the injection of urea, effluent geochemical data showed 
sustained production of NH 4 + . However effluent geochem- 
ical data also indicated that only a fraction of the urea 
added was hydrolyzed (-13%). Calculations based on the 
observed changes in Ca 2+ concentrations, when excluding 
effects of ion exchange, and alkalinity indicated that only 
0.12 millimole of calcium carbonate was precipitated per 
liter of solution used, which would have caused a negligi- 
ble increase in the surface area of the sediments packed in 
the column. Additional contribution to calcium carbonate 
precipitation from Ca 2+ exchanged from the sediment is 
small even when the exchange rate was assumed at the 
peak level observed (~ Day 16) throughout the experi- 
ment. Unaltered hydraulic conductivity measured at the 
end of the experiment supported the small amount of cal- 
cium carbonate precipitation with minimal impacts on 
hydraulic properties. 

The changes of & coincided linearly with the changes 
of g w throughout the experiment (Figure 4A). This linear 
correlation between & and o w agrees well with Archie's 
law [18] and indicated the dominant control of the pore 
fluid conductivity on the bulk electrical conductivity of 
the material. The most significant changes of & occurred 
at the beginning of the column experiment and the 



initiation of urea injection due to the significant changes 
in TDS associated with salt dissolution and/or ion deso- 
rption and urea hydrolysis, respectively. The slow 
decrease of both & and o w during the latter portion of 
the urea injection period coincided with the slow 
decrease of Ca 2+ , Sr 2+ and alkalinity (Figure 3A, B and 
3E), presumably due to the process of calcium carbonate 
precipitation. 

The observed changes of a" are more intriguing than 
the changes of &. As stated earlier, changes in the polari- 
zation (a") were expected during calcium carbonate pre- 
cipitation because of the additional contribution to the 
total surface area from newly precipitated calcium carbo- 
nate [29]. However, an increase of a" was not observed 
during the steady state of calcium carbonate precipitation 
in phase III, likely because the amount of calcium carbo- 
nate precipitation was small. This is in contrast to pre- 
vious experiments where a much larger amount of 
calcium carbonate was produced and strong electrical 
responses were observed [29]. It is important to note 
however that in the previous experiments large glass 
beads were used rather than natural sediments; the latter 
have a large baseline surface area that provides a strong 
and dictating background electrical polarization signal 
and further reduces the sensitivity of electrical signals to 
small amounts of calcium carbonate precipitation. 

Although not observed during the stable phase of cal- 
cium carbonate precipitation, we observed more signifi- 
cant changes in the polarization signal during the early 
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Figure 8 Modeled concentration profiles (lines) and measured data (squares) in the column effluent, for cases considering the full 
model (solid lines), suppressing ion exchange (dashed lines), and suppressing ureolysis (dotted lines). Time zero represents the end of 
Phase II at -13 days from the start of the experiment. 



stage of Phase I and III where significant changes of 
pore fluid chemistry occurred. Since o w changed signifi- 
cantly during these periods, it is reasonable to consider 
the changes of pore fluid chemistry as the possible 
cause of the observed changes in &'. 

With minimal changes to the sediment mineralogy, 
the observed changes in electrical phase response can be 
attributed to changes in EDL properties (e.g. charge 
density, charge mobility) in response to pore fluid che- 
mical alterations. This has been explored by previous 



researchers investigating the electrical properties of por- 
ous media [20,21,24,53]. We focus our discussion on the 
initial stages of phase I and III because the most signifi- 
cant changes in both polarization and fluid chemistry 
occurred during these periods. Figure 3 shows clearly 
the fluid chemistry differences during these two phases: 
while the major ions during Phase I were typical of site 
groundwater, new species were produced in Phase III, 
primarily NH 4 + , due to urea hydrolysis. Although urea 
was also a new species introduced during Phase III, we 
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Figure 9 Modeled calcium carbonate amounts within the column along its length (excluding the effluent) using (A) the full model and 
(B) without ureolysis. Amounts computed for the case without ion exchange differ only slightly from results of the full model (see Figure 8). 



assume minimal polarization effect from urea since it is 
presumably an uncharged species. During the early stage 
of Phase I, the concentrations of the major ions (for 
example Ca 2+ and S0 4 2 ~) decreased due to equilibration 
of the sediments with site groundwater. The initial high 
concentrations of the major ions were caused by salt 
dissolution and/or ion desorption from mineral surface 
as discussed above. Coincident with the decreases of the 
major ion concentrations was the continuous increase of 
the polarization signals {a"). At the beginning of phase 
III, the major pore fluid chemical changes were the pro- 
duction of NH 4 + and increase of alkalinity. As discussed 
above, the production of NH 4 + enhanced ion exchange 
and caused the occurrence of the short peaks for Ca 2+ 
and Sr 2+ . A concurrent increase of a" was observed dur- 
ing this period before reaching a steady state following 
the trend of change for NH 4 + . Comparison of the 
changes in both a" and fluid conductivity between 
phase I and III reveals an inversed correlation: the 
increases of a" in both phases, but the fluid conductivity 
was decreasing in phase I while increasing in phase III. 

The reversed correlation between a" and fluid conduc- 
tivity between the two phases suggests inherent differ- 
ences in the mechanisms of the polarization. The 
increase of a" is an indication of the increased capacitive 
effects at the grain/fluid boundary. Previous studies have 
indicated a binary change of this capacitive effect for dif- 
ferent ranges of ionic strengths; in particular, this capaci- 
tive effect was positively correlated with pore fluid ionic 
strength at low ionic strength, but negatively correlated 
at high ionic strengths [24]. The groundwater used in 
this experiments had an ionic strength of < 10 mM and 
generally falls into the low ionic strength range. Thus, 
the decrease of ionic strength (fluid conductivity) 
observed during the early stage of Phase I would have 
been predicted by Lesmes and Frye [24] to have resulted 



in a decrease of cr", contradictory to our observations. 
Insights learned from P wave data could potentially 
explain this discrepancy. The P wave velocities decreased 
systematically at the sensor levels sampling the native 
sediment during the early stage of Phase I. If a high frac- 
tion of swelling clays, specifically smectite, was present, 
the P wave velocity decrease could have been due to a 
gradual hydration and expansion of clay particles trapped 
at grain contacts. The existence of significant smectite 
fractions at the INL site has been previously documented 
[41]. Such a hydration process would provide additional 
polarizable surface area that was previously not available 
and thus promote polarization behavior, causing 
increases in &\ The fact that the zone of clean packing 
sand (level 4) did not exhibit the same velocity reduction 
suggests that a property change in the native sediment is 
responsible. The small but significant variation observed 
in baseline V p values at the different sensor depths was 
likely associated with heterogeneous packing of the sedi- 
ments within the column as mentioned above. 

The correlation between ionic strength (fluid conduc- 
tivity) and a" is positive during the early stage of Phase 
III, indicating a promoted polarization magnitude in this 
phase. We attribute this increase of c" during this stage 
to (1) production of NH 4 + that could potentially modify 
sediment surface charge properties through cation 
exchange and (2) hydroxyl ion (produced during urea 
hydrolysis) adsorption/titration on the mineral surface 
due to the pH buffering capacity of the sediments. 

As discussed previously, EDL polarization considers 
primarily the tangential movement of charges and the 
charge complexation structure of the EDL is expected to 
be relatively stable when pore fluid chemistry is unal- 
tered. However, active ion exchange behavior occurred in 
this system during Phase III due to the production of 
NH 4 + , an ion that is readily exchangeable with surface 
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complexed cations, including Na + , Ca 2+ and Sr 2+ . This 
ion exchange process could modify surface charge struc- 
ture, including charge density and mobility and could 
potentially result in an increase of the polarization at the 
interface. The effects of ion exchange on electrical signa- 
ture are less understood as has been observed and dis- 
cussed by others [22,23]. In addition, hydroxyl ions 
produced during urea hydrolysis have a tendency for sur- 
face adsorption/titration due to the high pH buffering 
capacity of sediments with large fractions of fine parti- 
cles. Hydroxyl adsorption/titration on sediment grains 
could enhance surface charge density and promote polar- 
ization behavior, as observed in other studies with porous 
media of high pH buffering capacity [54]. At the later 
stage of Phase III, ion exchange and hydroxyl adsorption/ 
titration presumably reached a steady state. This stabi- 
lized the surface charge structure with NH 4 + acting as 
one of the new ionic species contributing to surface com- 
plexation. Manifested in o'\ the polarization magnitude 
was stabilized during this period. 

The P-wave seismic measurements during phase III did 
not record an increase in velocity that would be indica- 
tive of grain cementation due to calcium carbonate preci- 
pitation. In light of the relatively small amount of CaC0 3 
precipitated and the localized nature of mineralization 
observed via micro- CT, this lack of a P-wave signature is 
not surprising. Previous observations of increases in 
frame shear stiffness during microbial urea hydrolysis 
[13] were associated with larger volume fractions of 
cement more evenly distributed across the pore space. 
The small changes in seismic response were likely driven 
by sample hydration and clay swelling rather than CaC0 3 
precipitation. 

The dynamic micro-CT observations of the micro-col- 
umn and the extracted sub-core clearly exposed evidence 
of mineralization. Since the precipitate observed was never 
directly interrogated using a method capable of resolving 
chemical composition, we have only indirect evidence that 
the mineralized patches were calcium carbonate; the com- 
bination of cation data and modeling suggests that this is 
likely the case. However, the occurrence of localized disso- 
lution/re-precipitation reactions involving silicate minerals 
within the sediment matrix cannot be ruled out, but 
would be expected to proceed at a much slower rate and 
thus to a much more limited extent than calcium carbo- 
nate, and would not be expected to significantly impact 
effluent chemistry over the relatively short duration of the 
test. The pore scale geometry of precipitates appears to be 
largely grain coating, localized in small patches rather than 
homogeneously distributed. This distribution, in addition 
to the minimal porosity reduction, likely contributes to the 
lack of observed changes in seismic P-wave velocity and 
hydraulic conductivity following urea hydrolysis. A large 
seismic signature would indicate a significant increase in 



mean contact stiffness, which would require precipitate 
formation at a larger number of grain contacts. 

After taking into account ion exchange and atmospheric 
equilibration processes, the latter of which induced prema- 
ture calcium carbonate precipitation and ammonia oxida- 
tion, the reactive transport model successfully simulated 
the concentration changes of major constituents during 
urea hydrolysis and calcium carbonate precipitation. The 
ion exchange processes associated with NH 4 + production 
were evidenced by effluent geochemical data and success- 
fully captured in the model. Changes in TDS calculated 
from variations of the concentrations (measured and pre- 
dicted) of major ionic species in the pore fluid correlated 
well with observed bulk electrical conductivity measure- 
ments; the change of both parameters was -10% - 12% 
during steady state of Phase III. Although the extent of 
urea hydrolysis was low (only -13% of the added urea was 
hydrolyzed), compared to predicted baseline levels the cal- 
cium carbonate precipitation was enhanced 3 to 4 times 
by the urea addition, confirming the successful promotion 
of calcium carbonate precipitation processes during the 
experiment. Research has indicated that attached micro- 
bial cells are responsible for a much greater fraction of 
urea hydrolysis than planktonic cells [8] . The low hydroly- 
sis rate in our column was probably associated with a low 
density of attached cells in the sediments based on cell 
count results presented above. As was reported earlier, the 
sediments used in the experiment were samples that had 
been stored dry in archive for over 10 years. Diminished 
microbial activity would have been expected. Although 
fresh groundwater was used to rehydrate the sediments, 
the attachment of cells to the sediments appeared to be 
limited based on cell adherence tests. Dissolution of Ag 
and AgCl from the electrodes used for geophysical mea- 
surements might have further inhibited urease activity [55] 
and have been observed by other researchers with extra- 
cellular urease (personal communication, George Redden, 
INL). 

Conclusion 

In this study, multiple geophysical methods including elec- 
trical, seismic and micro-CT, together with reactive trans- 
port modeling, were evaluated for their utility in 
monitoring and characterizing ureolysis promoted calcium 
carbonate precipitation for trace metal sequestration. 
Urease activity was successfully stimulated in columns of 
natural sediment and groundwater following the introduc- 
tion of molasses as a nutrient source. The rate of urea 
hydrolysis was low, probably due to the low density of 
microbes in the sediments and perhaps the toxicity of Ag 
from electrodes. However, this low ureolytic activity still 
led to a calcium carbonate precipitation rate that was 3-4 
times higher than the predicted baseline, supporting the 
potential of remediation approaches based on 
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ureolytically-driven calcium carbonate precipitation for 
sequestration of divalent trace metal and radionuclides. 
Although seismic monitoring was not able to detect 
changes during the calcium carbonate precipitation phase 
of the experiment due to minimal changes to the bulk 
modulus, complex conductivity measurements successfully 
detected changes in pore fluid chemistry during urea 
hydrolysis. Electrical conductivity correlated well with 
TDS changes due to urea hydrolysis and was successfully 
captured by effluent geochemical data and simulated by 
reactive transport modeling. Important for field monitor- 
ing of urea hydrolysis, we believe that the imaginary com- 
ponent of the electrical geophysical signal was related to 
the alteration of surface charge structure due to ion 
exchange from the production of NH 4 + and possibly 
hydroxyl adsorption to sediment surfaces or titration of 
surface proton sites. Previous observations of complex 
conductivity signatures directly from calcium carbonate 
precipitation [29] were not observed in this experiment 
presumably due to the small quantity of precipitated cal- 
cium carbonate, thus the small increase in surface area, 
relative to the volume of background sediments. Micro 
CT imaging demonstrated the sparseness of calcium car- 
bonate precipitation in the column and revealed the pat- 
chy and pore filling natures of the mineral precipitates. 
Micro-CT was shown to be a useful tool for dynamic ima- 
ging of the column to acquire direct evidence of precipita- 
tion geometry and distribution within the porous media. 
One dimensional reactive transport modeling with 
TOUGHREACT successfully simulated changes of major 
ion concentrations during the experiment, and the calcu- 
lated TDS changes of pore fluid based on the model and 
effluent geochemical data correlated well with electrical 
conductivity measurements. The reactive transport model- 
ing validated the reaction network outlined in Eq. 1-4 and 
demonstrated the usefulness of this model for simulating 
urea hydrolysis and calcium carbonate precipitation under 
field relevant conditions. 

Our experiments uniquely combined complex electri- 
cal, seismic, and micro-CT imaging technologies with 
reactive transport modeling to study urea hydrolysis 
coupled with calcium carbonate precipitation under field 
relevant conditions. Indirect geophysical measurements 
provided a quick indicator of system transitions between 
different phases while geochemical and modeling results 
offered detailed information on critical geochemical pro- 
cesses and parameters, such as ion exchange and 
changes in TDS. In addition to elucidating the evolution 
of the pore fluid geochemistry, the modeling was parti- 
cularly helpful in understanding intriguing geophysical 
signatures (e.g. imaginary conductivity) that were other- 
wise difficult to interpret. In addition, the precipitation 
end products inferred from all of these datasets were 
validated through micro-CT imaging that provided 



direct evidences of physical changes, including precipi- 
tate distribution and geometry. Although the limited 
extent of urea hydrolysis and calcium carbonate precipi- 
tation limited the demonstration of the complementary 
nature of these tools, our study showed clearly the 
advantages of combining multiple approaches to under- 
stand complex biogeochemical processes in the 
subsurface. 
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